arXiv: 1504.07605vl [astro-ph.HE] 28 Apr 2015 


Draft version April 29, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


THE FATE OF THE COMPACT REMNANT IN NEUTRON STAR MERGERS 

Chris L. Fryer 1,2 , Krzysztoff Belczynski 3 , Enrico Ramirez-Ruiz 4 , Stephan Rosswog 5 , Gang Shen 6 , Andrew W. 

Steiner'' 8 

Draft version April 29, 2015 

ABSTRACT 

Neutron star (binary neutron star and neutron star - black hole) mergers are believed to produce 
short-duration gamma-ray bursts. They are also believed to be the dominant source of gravitational 
waves to be detected by the advanced LIGO and the dominant source of the heavy r-process elements 
in the universe. Whether or not these mergers produce short-duration GRBs depends sensitively on 
the fate of the core of the remnant (whether, and how quickly, it forms a black hole). In this paper, 
we combine the results of merger calculations and equation of state studies to determine the fate of 
the cores of neutron star mergers. Using population studies, we can determine the distribution of 
these fates to compare to observations. We find that black hole cores form quickly only for equations 
of state that predict maximum non-rotating neutron star masses below 2.3-2.4 solar masses. If quick 
black hole formation is essential in producing gamma-ray bursts, LIGO observed rates compared to 
GRB rates could be used to constrain the equation of state for dense nuclear matter. 

Subject headings: Supernovae: General 


1. INTRODUCTION 


_Sin ce t heir discovery in the late 1970s (Klebesadel and 

Strong|1976 ), scientists have proposed a growing number 


of progenitors and engines for gamma-ray bursts (GRBs). 
The variety in the observations of classical GRBs points 
to multiple progenitor scenarios. Based on their dura- 


GRBs ( 

Kouveliotou et al. 

1993) 

accretion disk (BHAD) c 

lass o 


1993). Within the black hole 


bursts are believed to arise from accreting black h oles at 


the ce nter of massive stars, e.g. from co llapsars (Woosley 


19931 or he-mergers (Fryer et al. 19991. Short-duration 


bursts are believed to be produced by the merger of two 
compact objects: consisting either of a binary neutron 
star (NS-NS) or a neutron star and a black hol e (BH - 


NS): for review see, f o r example, Fryer et al. 


Popham et al. (19991; Lee and Ramirez-Ruiz 


(1999); 

(2007); 


Rosswog (|2U15[) . These compact mergers are produced in 


close binaries and the kicks imparted onto these compact 
remnants cause these systems to have high space veloci¬ 
ties. These velocities, coupled with merger times, mean 
that bursts produced from NS-NS and BH-NS binaries 
should have broad spatial distributions with respect to 
their host galaxy, including bursts that occur well out- 
side of their host (Bloom et al. 1999 Fryer et al. 1999 
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Zemp et al. 2009). With Swift, the number of localized 
short-duraUon Bursts has increased dramatically, and the 
observed spatial distribu tions of these bursts matches the 


predictions from theory (Belczynski et al. 2006 

Fong and 

Berger 

2013 

Bchroozi et al. 

2onp indiBa; 

no other 


model to date can easily explain the spatial distributions 
of short-duration bursts, and NS-NS and BH-NS mergers 
are almost universally considered the leading progenitors 
for these bursts. 

The BHAD engine is not the only way accreting 
compact objects can produce possible outflows. Neu¬ 
tron s tars are also able to accrete at Super-Edd i ngton 
rates (Houck and Chevalier 1991 Fryer et al. 19961 
and NS accretion disk (IN SAL)) systems will look very 
similar to BHAD engines. Alternatively, building off 
of the leading mechanism for soft gamma-ray repeaters 
(SGR), theori sts have a rg ued t hat rapidly-spinning mag- 
netars dDun can and Th ompson 1992 )_can produce clas¬ 
sical GRBs (|Zhang and iVLeszaros 2001). The disadvan¬ 
tage of the neutron star models is that they - as long 
as the merger remnant is stable - can driv e stro ng bary- 


onic winds via neutrino energy deposition (Dessart et al. 
20091 Perego et al |2014|) that potentially choke the GRB 
Ramir 


let (IRosswog and Ramirez-Ruiz||2002[ |Murguia-Berthier| 
et~al.||2014| ). 


The advantage of the neutron star model is that the NS 
can be born with strong magnetic fieldt|^] and these mag- 


1992 

Thompson 

1994 

Rowlinson et al.j 

2014 

). Under- 

standing the relative rate ot merging sys 

bems 

that form 


black holes versus neutron stars can help direct theorists 
toward a better understanding of these systems. For ex¬ 
ample, if very few systems collapse to form black holes, 
either the baryonic loading problem is less severe than 


9 Note that generally in high accretion scenarios like these merg er 
systems, there is a belief that the accretio n buries the field(|Popov| 
and Turolla 2012 Vigano and Pons 2012). Typical timescales lor 
the re-emergence ot these magnetic fields are tens of kyr. In such 
a case, any magnetar-like model will not work to explain a GRB. 
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currently believed or mergers are not the solution to 
GRBs. On the other extreme, if only a small fraction 
remain neutron stars, it is worth identifying and trying 
to observe the failed GRBs formed by these systems. Ob¬ 
servations of relative rates of these systems can provide 
a fairly direct means to constrain the nuclear equation of 
state. 

This paper brings together studies of binary neutron 
star merger models, nuclear equations of state, and pop¬ 
ulation synthesis models to study the formation rate 
of different potential GRB progenitors, differentiating 
between black hole accretion disk models, accretion- 
induced collapse models, and neutron star accretion 
models. Section [2] describes the merger models used in 
this study and Section [3] describes the effects and un¬ 
certainties in the equation of state (EOS). By coupling 
these results with population synthesis studies, we can 
compare these results to the suite of GRB observations 
(Section®. We conclude with a discussion of the impli¬ 
cations of these results on the wide range of phenomena 
explained by NS-NS and BH-NS binaries. 


2. MERGER MODELS 

In the past decade, the number of hydrodynamic simu¬ 
lations of NS-NS and BH-NS mergers has grown consid¬ 


erably (|FbucarFeFaI!p010|[20TT||Foucartp012| Korobkin 

.|2(jl2| I'Hokokczaka et al.|2(Vl:]l |Ryut 0 ku 


etlir 2013; 


et al. 


Takami et al. 2014; Bauswein et al. 2014 

IRadice et al. 

2014 

Shibata et a.|2014 Kiuchi et aJ2014 

Foucart et al. 


sophisticated, most of the current models include only a 
subset of the physics needed to model these objects: hy¬ 
drodynamics, neutrino transport (or at least neutrino en¬ 
ergy losses), nuclear equations of state, magnetic fields, 
and general relativity. Nevertheless, these models are 
gradually painting a complete picture of the merger pro¬ 
cess and we can use them to make a first pass at the fate 
of these systems. 


2.1. BH-NS mergers 

To date, no BH-NS system has been observed and, at 
this point, there are no observations that would require 
such systems to exist. However, there are clear observa¬ 
tional biases against such systems: if the NS is formed 
after the black hole, as we expect, it is difficult to re¬ 
cycle it. Hence we do not expect a millisecond pulsar 
in BH-NS binaries. Population synthesis models predict 
that these systems form frequently and many of them 
close enough to merge within a Hubble time. Theoreti¬ 
cal models for massive star evolution and the formation 
of BH-NS systems suffer from large uncertainties (e.g., 
treatment of star internal mixing, stellar winds, common 
envelope and supernova explosion). Population synthe¬ 
sis models predict a range of formation rates from several 
orders of magnitude below the BN S me rger ra te to rates 
that rival the BNS rate (Dominik et al.||2012|. 

The binary parameters of BH-NS systems are crucial 
for the outcome of the merger and only for a fraction 
of the parameter space will it be possible to form an 
accretion disk massive enough to launch a GRB. If we 
assume the energy arises from the disk, 


E P: 


2 x 10 51 erg ( — 


(c) 


-Mdisk 

O.OIMq 


( 1 ) 


can be extrated, so that the disks do not need to be 
extremely massive to accomodate the typ ical isotro pic 
gamma-ray energies of short bursts (Berger 2014) of 
~ lO 50 erg (and correspondingly lower if they are colli¬ 
mated) . If the energy arises from the rotating black hole, 
significantly smaller disk masses are needed to pr oduce 


the energie s required for short bursts (Lee and Ramirez- 
Ruiz||2007 ). Nevertheless, not all systems will be able to 


form disks, since the radius where tidal disruption sets 
in, must be larger than the innermost stable cir¬ 

cular orbit (ISCO), -Risco- Since Risco °c Mbh, but 

Rtid oc M b ' h , for massive enough black holes the tidal 
radius lies inside of the ISCO and forming a massive 
accretion disk becomes impossible. For a non-spinning 
black hole this occurs already near Mrh = 8 Mm, so 
that BHs of the masses th at are thought (Belczynski et| 
al. 2008 Ozel et al. 2010) to be most likely (~ 10 M 0 ) 
need large dimensio nless spin para meters, aBH > 0.9, to 
form sizeable disks (|Foucart 2012). Finally, the inclina¬ 
tion of vector of BH-JNS orbital angular momentum to 
BH spin vector plays an important role in the forma¬ 
tion of a torus. At large inclination angles (> 40 — 90°), 
material from the disrupted NS is ejected from the vicin¬ 
ity of a merger and does not form torus. Non-zero tilts 
are expected in BH-NS systems if a NS is formed as a 
second compact object in a binary and with a non-zero 
natal kick in isolated binary evolution ffifield popula¬ 


tions (Fryer et al. 1999 Dominik et al.|20~12 2013 2014) 


globular clusters! 

_--2 

Grindlay et al. 

2006) l.ee et al. 2011): 

East et al. 

2013 Rosswog 2013; 

Tsang 2013; Ramirez- 

Ruiz et al. 

2015). 


their fates has grown with time ( 

Rosswog 2004 2005[ 

Rantsiou et al. 2008; 

Foucart et al. 

2010, 2011; Kyutoku 

et ai.|2013 Paschaiid 

is et al. 2014[ 

Foucart et al. |2015|), 

all showing that the 

tinal fate ot tr 

le merged system de- 


pends sensitively upon the initia l conditions. Based on 


earlier merging BH-NS models by Rantsiou et al. (2008) 
it was shown for one particular evolutionary model that 
only a fraction of BH-NS systems may potentially form 
an accretion torus and lead to a GRB. For small BH spins 
(a sp in < 0.6), only 1% of BH-NS mergers were found to 
have a disk needed to produce a GRB. For high BH spins 
the fraction o f systems with di s ks hav e been calculated 
to be ~ 40% Belczynski et al.| (|2008|). It is important, 
that in future studies, the merger hydrodynamical simu¬ 
lations are mapped with astrophysically-motivated pre¬ 
dictions of BH-NS system physical properties. 

For this paper, we will assume the GRB rate is dom¬ 
inated by BNS mergers, but we will return to the topic 
of BH-NS mergers when we better understand the BNS 
results. 


2.2. NS-NS mergers 

NS-NS mergers always form sizable accretion disks and 
aare the more canonical GRB model. The primary diffi¬ 
culty with these mergers as GRB progenitors is the na¬ 
ture of the merged core. After the merger, the remnant 
will evolve according to four separate pathways: collapse 
directly to black holes, those that initially form NSs but 
subsequently collapse during disk accretion, those that 
don’t collapse to black holes until after the disk has fully 
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accreted and the newly formed neutron star spins down, 
and those that, even after the spin down, remain a NS. To 
determine the fate of each merged object, we must com¬ 
bine both simulations of the merger process (discussed in 
this section) with our current understanding of the nu¬ 
clear equation of state (see Section |3|. For two equally- 
massed neutron stars, the merger produces a condensed 
core consisting of two major components formed from 
the cores of the initial neutron stars. As the mass ratio 
becomes more extreme, the more massive neutron star 
remains more intact, disrupting its companion. In all 
cases, the merger produces a merged core surrounded by 
a dense accretion disk with a small amount of material 
(~1%) ejected during the tidal ejection process. 

To study the fates of neutron stars, we use the large 
grid of smooth particle hydrodynamics simulations from 


Korobkin et al. (2012). This grid of models includes a 
broacTsuite of neutron star mass pairs using the same in¬ 
put physics. For each simulation, we calculate the mass 
and angular momentum of the merged core (all mate¬ 
rial above 10 14 gcm“ 3 ), disk masses (bound mass with 
densities below 10 14 gcm -3 ) and ejecta masses (TableJT]) . 

At the end of the smooth particle hydrodynamics sim¬ 
ulations, considerable mass (~ 0.5 — 0.9 Mq) remains 
bound but at densities below 10 14 gcm -3 . In Table[l] we 
refer to this material as “disk” material and much of it 
will form a disk that quickly accretes onto the compact 
core. If the core collapses to a black hole, the system 
immediately evolves into a BHAD engine. If not, it may 
first pass through a neutron star accretion disk or mag- 
netar engine until the accretion drives the neutron star 
above the maximum mass and it collapses, producing a 
BHAD scenario. 

Whether or not the core is a black hole or neutron 
star depends upon whether it is more massive than the 
maximum neutron star mass at its spin rate. This de¬ 
pends upon the still-uncertain equation of st ate. Recent 
observations by Demorest et al. (2010) and Antoniadis 
et al. (2013) demonstrate that this mass is least 2M 0 . 
The exact value for the maximum mass of a particular 
merged core depends on the equation of state as well as 
the internal energy and angular momentum distribution 
in the core. Our choice of equation of state will deter¬ 
mine which merged cores collapse immediately to a black 
hole. 

In other cases, the core must accrete some material 
before it collapses. To determine the engine behind 
this neutron star, we must understand the accretion 
timescale. If the core remains a neutron star for a long 
period of time, the neutrino-driven wind will choke the 


ou tflow, preventing a strong burst (Murguia-Berthier et 
aL|2014). To estimate the accretion time of the disk, we 


use the matter distribution from our models and assume 
an a-driven disk scenario. That is, for the spatial dis¬ 
tribution of the disk material, we assume the accretion 
timescale is just the orbital period divided by the value 

« 3/2 

of a: T acc = ^ (gm qi /2 where r is the spatial position 
of the matter, M enc \ is the enclosed mass and a is the 
viscous disk parameter. Figure [I] shows the mass of the 
core as a function of time under our a-disk assumptions 
using a = 0.01. Note that although much of the accre¬ 
tion occurs in the first 100 ms, the accretion phase can 
last out to a few seconds for some systems. Clearly, the 


accretion time depends linearly on the value for a whose 
value is uncertain to an order of magnitude. 



Log Time (ms) 


Figure 1 . Core mass versus time for our merged systems assum¬ 
ing an a-disk solution for the accretion time. Note that in many 
systems, the accretion occurs very rapidly (in the first 10,ms) and 
can last out to a few seconds. The jumps occur because the an¬ 
gular momentum of material is not uniform and large amounts 
of material have similar angular momenta, and hence, accretion 
timescales. 


rl|1992 Kochanek||1992 ). 
ibstantially smaller than 


For the most part, these models assume that the neu¬ 
tron stars in the binary are initially not spinning (al¬ 
though for one case, we include a neutron star binaries 
in co-rotation). The assumption is well-justified since 
the last inspiral stages happen so quickly that even an 
implausibly l arge dissipatio n cannot spin u p the star s 
substantially (Bildsten & Cutl er|| 

Pre-merger spins are typically substantially i 
co-rotation, but even in the co-rotation case, the mass is 
altered by only a few percent, but the angular momen¬ 
tum in the core can change by 15%. This has only a 
small affect on the maximum stable neutron star mass 
of the core, but it can dramatically change the accretion 
rate (see Section |3h . 

Be aware that tire current set of simulations are in¬ 
complete, missing both physics and spatial resolution to 
model the physics. As this physics is improved, exact 
quantities on merged cores, etc. will move. But the ba¬ 
sic trends will continue and, although there may be a 
shift in our results, the trends and methodology to study 
NS-NS and BH-NS mergers the results will remain. 


3. EQUATION OF STATE 
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The fate of the core depends on the exact details of 
the equation of state. If the maximum neutron star 
mass were 2.0M f T 1 , the maximum observed neutro n star 
mass (Demorest et al.|2010 Antoniadis et al.|2013 ), most 
of our merger cases would immediately collapse to a black 
hole. This observation only places a lower limit on the 
maximum mass, and the maximum mass of a cold neu¬ 
tron star may, depending upon the equation of state, be 
much higher than this maximum observed mass. In ad¬ 
dition, fast-rotation and thermal energy can significantly 
raise the maximum mass. To determine the final fate of 
our merger remnants, we will study the maximum mass, 
and its dependence on the spin and internal energy, for 
a range of equation of state models. 

The merged core, especially if the components are 
near-equal mass, are often born with strong differential 
rotation (Rosswog and Davies|[2002). Depending upon 
the magnetic held strength, these cores will v ery quickly 
evo lve in to uniform rotation. For example, |Shibata et| 
faT) (2007) argue that a combination of magnetic wind¬ 
ing, shearing and the magnetorotational instability can 
quickly redistribute the angular momentum in an accret¬ 
ing disk system. For a few of our fastest rotating mergers 
(the merger of a 2.0 M 0 NS with a 1.0,1.2 M 0 NS and the 
merger of a 1.8 M 0 NS with a 1.0 M 0 NS: see Table [l]), 
secular instabilities can develop that will also redistribute 
the angular momentum. Finally, neutrin o cooli ng in 


th e core will drive convection in the core (Rosswog et 


al. 2003) that can redistribute the angular momentum. 
These redistribution mechanisms have timescales lying 
between 10-100 ms. If it is at the low end, our assump¬ 
tion that it is nearly instantaneous holds. If there is a 
longer delay in this evolution to uniform rotation, the 
maximum mass could remain high for longer, produc¬ 
ing more long-lived neutron star systems. Note, how¬ 
ever, that for our population study (see section EL only 
0.3% of our systems have a NS component with masses 
above 1.8 M 0 and none of our ~ 6000 systems included 
pairs with a massive component above 1.8 M 0 and a low 
mass component below 1.2 M 0 , so the highest angular 
momentum systems are extremely rare. 

Although a multitude of dense matter equations of 
state (EOS) exist which make firm p redictions abou t 
maximum masses, see for example, Cook et al. (1994), 
all of them make simplifying assumptions for the physics 
and theory can not place strong constraints on the fi¬ 
nal maximum mass. Most models for the neutron star 
equation of state predict masses in a range between 2.0 
and 3.OM 0 , uncertainties in nuclear physics prevent a 
more precise prediction of the maximum mass for a non¬ 
rotating model. However, for a model with a given max¬ 
imum non-rotating mass, we can predict the effect of 
thermal energy and angular momentum. 

For our study, we use a variety of equations of state 
to probe the fate of our merged cores. Thes e in¬ 
clude, in particular, the NL3 (Lalazissi s et a l. |1997l and 
FSU2.1 (Todd-Rutel and Piekarewicz||2005p EOSs. We 
also in clude one of the equation of sta te models devel¬ 
oped in 
pararne 


Steiner et al. (2010 2013 2015) where the model 
ers are matched to observations of quiescent low- 
mass X -ray binaries and photospheric radius expansion 
bursts (Steiner et al.||2015) usin g a Monte Carlo schem e 
to perform marginal estimation ( Steiner et al. |2014b|c ). 
In this equation of state model, matter near and below 


the nuclear saturation density is described by a parame¬ 
terization calibrated to the properties of labor atory nu¬ 
clei and quant um Monte Carlo simulations (Steiner and 
Gandolh 2012). At higher densities, matter is described 
by two piecewise polytr opes. Th e RNS code for rotating 
neutron stars from Stergioulas and Friedman (1995) was 
embedded in our Monte Carlo simulation to output prob¬ 
ability distributions for the maximum mass of a neutron 
star given a fixed angular momentum. 

During the merger, shocks and viscous forces increase 
the energy in the core, raising the entropy. This ex¬ 
tra th ermal energy can ra i se the maximum neutr on star 
mass (Kaplan et al. 2014 Bauswein et al.| 2014). Most 
studies where this has been a large effect have assumed 
entropies of roughly 8 fee per nucleon. However, ac¬ 
tual es timates of the post-me rger entropies used in our 
study flKorobkin et al.||2012|) lie at about 1 kn per nu- 
cleon (Table flj) . Using the XL.8 (jbalazissis et al.| |1997 I 
equation of state (fairly typical) an entropy of 1.28 fee 
per nucleon raises the maximum stable mass by only 1%. 
In this scenario, the internal energy deposited in the core 
in the merged system is not sufficient to strongly affect 
the maximum neutron star mass. There is more ther¬ 
mal energy in the surrounding disk, but its low densities 
allow efficient neutrino cooling that is on par with the 
accretion time. Hence, we will use the accretion time to 
estimate our collapse rate for systems that collapse after 
accretion. 

The angular momentum in the merged core can, how¬ 
ever, have a larger impact on the maximum neutron star 
mass. If we assume the core quickly evolves into a sys¬ 
tem in solid body rotation, we can use existing equations 
of state to estimate the maximum neutron star mass for 
these rotating systems. For the FSU2.1 equation of state, 
we find that the maximum neutron star mass increases 
nearly linearly with rotational support over the range of 
rotations produced in neutron star mergers (Fig. |2]). As 
input for our population studies, we have fit the max¬ 
imum neutron star mass predicted for this equation of 
state with the following linear approximation: 

M max , NS = (2.1 + J/(10 50 g cm 2 s _1 ))M 0 (2) 

where J is the rotational angular momentum of the com¬ 
pact core. Figure [2] also shows the distortion of the neu¬ 
tron star that would form with such high angular mo¬ 
menta and this distortion reflects the large asymmetries 
produced in the cores of these merger events. 

Thus far we have only used specific equations of state. 
To determine the role of angular momentum change with 
a broad range of allowed nuclear equations of state, we 
use our parameterized equation of state. The increase 
in the mass (AM) of a maximally-rotating neutron star 
depends upon the maximum mass of the non-rotating 
neutron star. Figure [3] shows the distribution of AM 
as a function of maximum non-rotating mass. For a 
2.0 Mq star, a maximally-rotating neutron star lies be¬ 
tween 2.35 — 2.5 Mq depending upon the exact equation 
of state. For an equation of state producing a maximum 
non-rotating mass of 2.6 M 0 , this mass is 3.16 — 3.22 M 0 , 
an increase of 0.56 — 0.62 Mq. The AM for a maximally 
rotating neutron star increases for equations of state that 
have larger maximum non-rotating stars. However, for 
the angular momenta in our merged cores, the increase in 
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Figure 2. Maximum neutron star mass versus angular momentum 
for the NL3 equation of state for a set of central pressures (top). 
A neutron star with a maximum non-rotating mass of 2.1 Mq will 
have a maximum mass of 2.6Mq. This angular momentum can 
drastically distort the neutron star. The ratio of the polar to 
equatorial radii (assuming solid-body rotation) is also shown as 
a function of angular momentum. At the end of our simulations, 
the core structure has not yet reached solid body rotation, and 
the asphericity is even more extreme (possibly arguing for a more 
extreme maximum mass). 

the maximum mass is much less dramatic. A fit to these 
lower-angular momentum systems yields the following re¬ 
sult for the maximum rotating neutron star mass: 

M maXjNS = M 0 + Ml x (J/(10 49 gcm 2 s _1 )) /3 M 0 (3) 


where Mq is the maximum mass for a non-rotating 
neutron star. For this paper, we study different val¬ 
ues of Mo, fitting the change in mass to average val¬ 
ues of rotating systems for a range of angular momenta 
from 0 — 3 x 10 49 gcm 2 s _1 based on the angular mo 


'IT 

ia £ 


For 
f3 are: 


menta in the cores in our simulations (Table I 
Mq = 2.0,2.5,2.7M 0 , the values for Mi and 
Mi = 0.0219,0.0354,0.0574 and /? = 1.75,1.5,1.6 respec¬ 
tively. Typically, rotation adds less than a few tenths of 
a solar mass to our maximun mass of our merged cores. 
Accretion of the disk material can also alter the angular 
momentum of the core, but how much angular momen¬ 
tum depends upon the exact accretion process (including 
magnetic field strength). For this study, we assume the 
specific angular momentum is held constant during the 
accretion phase. Moreover, we assume that the post¬ 
merger magnetic fields are small enough (F?|||10 17 G) to 
not contribute noticeably to the pressure. To obtain a 
more detailed survey of the possible merger fates, we 
created a grid of maximum masses, varying Mq from 2.0 
to 2.5 with constant values for Mi (= O.O3M 0 ) and /? 
(=1.6). 

With these fits to the maximum mass of our equation 
of state, we can determine the fate of the merged systems 
presented in Sec tion [2] Here we use the four fates dis¬ 
cussed in Section |2.2| BH, BH acc , BH sp i n , NS. Note that 
for these fates, our equation of state models assume that 
the neutron star has reached an equilibrium density state 
and it is in solid body rotation. The rotating object is 
differentially rotating through this period but our simpli¬ 
fied assumptions provide an estimate of the importance 
of spin in the final fate. Spin may prevent an immediate 
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Figure 3. Probability distribution of AM(Mq) as a function 
of the maximum non-rotating mass. The top panel focuses on 
masses near and above a maximum mass of 2.0 Mq whereas the 
bottom panel zooms in on those maximum non-rotating masses 
above 2.5 Mq. 


collapse to a black hole but, in most cases, the accretion 
of the surrounding material is sufficient to cause the core 
to collapse with or without rotation (see Table [2]). 


4. TYING IT ALL TOGETHER: POPULATION SYNTHESIS 

To compare these results to observations, we must in¬ 
corporate our findings into a study of the binary neutron 
star and neutron star/black hole populatio ns. We employ 
the StarTrack population synthesis code (jBelczynskiet 
al. 2002 2008) to generate a population of binary com¬ 


pact objects (in particular NS-NS and BH-NS binaries). 
The code is based on revised formulas from Hurley et 
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al. (2000); updated with new wind mass loss prescrip¬ 
tions, calibrated tidal interactions, physical estimation 
of donor’s binding energy (A) in common envelope calcu¬ 
lations and convection driven, neutrino enhanced super- 
no va engines. A full desc ription of these updates is given 


Dominik et al. (2012). The two most recent updates 


take into account measurements of initial parameter dis¬ 


tributions for massive O stars (Sana et al.||20121 as well 
as a correction of a technical bug that has limited the 
formation of BH-BH binaries for high metallicity (e.g., 
Z = 0.02). 

Two major factors shape the final compact object 
mass: wind mass loss and core collapse/supernova com¬ 
pact object form atio n. For w ind mass loss we use O/B 
type winds from Vink et al. (2001) and for other evolu¬ 
tio nary stages (e.g., -L BV w inds) formulae as calibrated 
in IBelczynski et al.l (|2010 |). We adopt set of models 
presented by |.Fryerewir (2012) with the rapid core 
collapse/supernova mechanism. The explosion occurs 
within the first 0.1 — 0.2s driven by a convection and neu¬ 
trino enhan ced engine. This “rapid ” supernova engine 
reproduces (Belczynski et al. 2012) the observed mass 
gap (apparent lack of neutron stars and black holes with 
mass in range 2-5 M^) in Gal actic X-ray binaries (Ozel 
et al-polo Bailyn et al.||1998 ). Additionally, we include 


INS star formation through electron capture supernovae 
(ECS) for low mass progenitors M zams ~ 7-11 M 0 . The 
range of initial progenitor mass on Zero Age Main Se¬ 
quence (ZAMS) for ECS depends sensitively on interac¬ 
tions (mass loss and gain) of stars in binary. We assume 
that all NSs formed through ECS have mass of 1.26Mg. 
Stars initially more massive than ~ 10 Mg typically form 
NS in Fe core collapse through our adopted rapid super¬ 
nova engine. NS are formed in a broad mass range 1.1— 
2.5 Mg, although most NSs that form close NS-NS bina¬ 
ries are found in much narrower range 1.2-1.4 Mg with a 
long tail extending to 2 Mg . In our evolutionary calcula¬ 
tions we assume that ECS NSs do not receive natal kicks 


(Podsiadlowski et al. 2004), while for NS formed in Fe 
core collapse we adopt a 1-D Maxwellian with er = 265 
km s” 1 (t hat corresp onds to an average kick of ~ 400 
km s -1 (Hobbs et al' 2005). 

To generate the populations of double compact objects 
we adopt power-law initial mass function (IMF) with 
slope of —2.7 for the massive primaries (M zams = 5- 
150 M 0 ). Secondary mass (less massive binary compo¬ 
nent) is drawn from flat mass ratio distribution. Then 
we adopt two very different models for initial orbital sep¬ 
aration and eccentricity. In one set of models we adopt 
flat in logarithm (oc 1/a) distribution of initial separa¬ 
tions with thermal eccentricity distribution (= 2e). Such 
assumptions have been adopted in many population syn¬ 


thesis studies (e.g., for references and examples see Bel¬ 
czynski et al. (2008)). We refer to the models employing 


these assumptions as “OLD”. It was recently pointed out 
that these distributions may not be adequate for massive 
progenitors of NSs and BHs. We therefore adopt revised 
distributions: closer orbital peri ods oc l og(P)~°' 55 and 
less eccentric orbits oc e -0 42 ( {Sana et a l. 2012). We refer 
to models with these initial conditions as “1NEW”. 

Each initial condition calculation is performed for two 
characteristic stellar metallicities: solar composition Z = 
Z Q = 0.02 and Z = 0.1 Z Q = 0.002. We refer to these 


models as high and low-metallicity models. In our map¬ 
ping of population synthesis calculations with NS-NS 
merger calculations, we use evenly mixed (50% -50%) 
population of low and high metallicity models. This is 
to mimic (in very approximate way) a stellar content of 
Universe at various redshifts. 

In each calculation we allow for a different outcome of 
Roche lobe overflow (RLOF) initiated by a Hertzsprung 
gap (HG) donor. If RLOF appears to proceed on a dy¬ 
namical timescale we perform CE envelope energy calcu¬ 
lation to check whether a given binary survives the event 
(model A). Under this prescription, we allow for the effi¬ 
cient formation of close double compact object binaries. 
For contrast, we assume that each event of dynamical 
RLOF with HG donor does not lead to the formation 
of close double compact object (model B). Some of HG 
donors (right after main sequence) may not have devel¬ 
oped clear core-e nvelope structure and are a subject to 


CE merger (e.g., Belczynski et al. (2007)) or at later 
times HG donors may not have yet developed convective 
envelope and instead of CE they may potentially initiate 
thermal-timescale mass transfer that does not allow for 
significant orbital contraction and close double compact 
object formation. 

Additionally, we allow for uncertainties in our calcu¬ 
lation of NS mass. In particular we allow NS mass to 
increase by a fixed factor and in particular we chose this 
factor to be 0.1 Mg. Predictions of neutron star masses 
fold in a number of uncertainties from understanding the 
mass of the core at bounce, the duration of the engine 
phase and the amount of fallback. At the low mass end, 
fallback is not a concern an d the duration of the engine 
must be short (Fryer et al. |2012 ). But this initial NS 
mass can be larger by as much as 0.1 Mg. In many cases, 
the first neutron star formed also goes through a com- 


mon envelope phase where it can 
through hypercritical accretion I 

accrete up to ~u.i 
Fryer et al. 1996| 

iV U\ 

Bel- 

czynski et al. 2002, 2010b Mac! 

Leoad & Ramirez-] 

duiz 


20l5:1|j. 


Physical properties of double compact objects of all 
sorts (NS-NS, BH-NS and BH-BH systems) generated 
under above conditions are discussed in detail in a series 


of rec e nt papers IDominik et al.| (|2012|), |Dominik et al. 
(20131, Dominik et al. 12014|) ana de Mink fe Belczynski 


(2015, to be submitted) . The actual population synthesis 
models are available at www.syntheticuniverse.org. 
In particular, generated NS-NS merger rates (~ 10-40 
Myr -1 ) are in agreement with the empirically esti mated 
Galactic NS-NS merger rates (~ 10-100 Myr” 1 ; |Kim et 


al. (2015)). The NS masses predicted by our models are 


in general agreement with the mass estimates available 
from close NS-NS systems ( Fig. PH ). H owever , note that 


the X-ray binary 4U1538-52 van Kerkwijk et al. (1995) 


suggests that neutron stars are born with mases below 
1.1 Mg, suggesting the minimum neutron star mass can 
not be much higher than our standard value, echoing the 
theoretical error bars of O.lMg. 

With these models, we can determine the fate dis¬ 
tributions of merging systems (if and when they form 
black holes) as a function of the equation of state. Ta¬ 
ble m shows the results for a a series of models: with 
and without systems that undergo a common envelope 
in the Hertzprung gap, with the standard initial condi- 
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NS mass [M Q ] 



Figure 4. The neutron star mass distribution (number per mass 
bin as a function of mass) for our models using the “New” ini¬ 
tial conditions. These results include systems that have under¬ 
gone common envelope evolution in the Hertzprung gap with an 
enhanced neutron star mass (increasing it by O.IMq). For com¬ 
parison, we plot the masses of the observed NS-NS systems. 


tions and the new (Sana et al. 2012) initial conditions, 
with our standard and modified (increased by O.IMq) 
neutron star masses. The general trends between all 
these models are the same. The primary fates of the 
merged binaries are bimodal: either the system collapses 
to form a black hole after some disk accretion or the core 
remains a neutron star. Most of the systems that col¬ 
lapse to black holes do so within the first 100 ms. This 
means that most of the black-hole forming systems could 
also produce GRBs (avoiding the baryon contamination 
problem). 

The fraction of systems that collapse to a black hole 
is above 90% for any EOS that has a maximum non¬ 
rotating mass below ~2.3-2.4Mq, but drops precipi¬ 
tously beyond this and most models predict that over 
90% of the cores remain neutron stars for any EOS with 
a higher maximum rotating neutron star. This strong 
sensitivity to the EOS means that if we can observe these 
differences, we have a strong probe of the neutron star 
equation of state. We will discuss the observational im¬ 
plications of these fits in Section [5j 

Except at the boundary between most systems form¬ 
ing black holes and most systems forming neutron stars, 
the differences between the primary fates of our popu¬ 
lation syntehsis models are fairly small (few %). At the 
boundary, differences in the population of 10% are pos¬ 
sible. The primary uncertainty in our calculation is the 


lower limit on the neutron star formation mass. For our 
lower NS formation limit of 1.1 M 0) the Steiner2.4 (max¬ 
imum non-rotating neutron star mass of 2.4 Mq) marks 
the transition between mostly BH final states vs. mostly 
NS final states. If the mass is increased by O.IMq, the 
Steiner2.5 marks this transition. More importantly for 
our analysis is the transition EOS for systems that form 
BHs in less than 100 ms. For these systems, the tran¬ 
sition is Steiner2.2, Steiner2.3 respectively for our two 
mass limits. The other population synthesis parameters 
do not produce much variation in these results. 

5. IMPLICATIONS 

In our analysis, we found that the fate of the merged 
system depends sensitively on the neutron star equa¬ 
tion of state. Figure [5] shows the fraction of rapid 
(f acc < 100 ms) and immediate collapse black holes as a 
function of the equation of state. If we can observation- 
ally distinguish between these fates, we can use NS-NS 
mergers to place constraints on the equation of state. 



Figure 5. Fraction of systems that produce standard BHAD 
GRBs and Neutron stars for our 4 basic population models: stan¬ 
dard NS masses (BHAD - solid, NS - dotted) for old (blue) and new 
(cyan) initial orbital parameters, NS masses increased by O.IMq 
(BHAD - dashed, NS - long-dashed) for old (red) and new (ma¬ 
genta) initial orbital conditions. The orbital conditions have vir¬ 
tually no effect on these results and the variations caused by other 
binary population synthesis parameters is also negligible. The pri¬ 
mary error in these estimates is that of the initial mass distribution 
of neutron stars, producing an equivalent error in the maximum 
non-rotating NS mass limit on the EOS. 


For example, the current favored model for short- 
duration GRBs is the merger of two neutron stars form¬ 
ing a black hole accretion disk. If the merged core forms a 
neutron star, neutr ino-driven winds ( Dessart et al.|2 009 


XX\_. UU1 WXX O c CXI . XXV LX Ui XXXU UilVLXX VV XXXVXO V CXX. 

Perego et al.||2014| drive outflows through which any jet 


must penetrate to produce a gamma-ray burst. Current 
studies argue that this well-known baryon contamination 
problem can only be avo ided if the ne ut ron star phase 
lasts less than ~ 100 ms (Murguia-Berthier et al.||2014). 
















































If this is true, GRBs are produced in only a fraction 
of all merging systems. If the maximum mass is below 
2.3-2.4 M 0 , then the cores of most neutron star mergers 
produce black holes within ~ 100 ms of the merger. For 
these equations of state, most neutron star mergers are 
capable of forming short-duration GRBs. 

Above a maximum non-rotating neutron star mass of 
2.3-2.4 M 0 , less than 4% of all mergers produce rapid- 
collapse black hole systems within 100 ms. If only such 
black hole systems form short-duration GRBs, this would 
require a rate of mergers that is 50 times higher than 
cannonical values. Although possible, this strains our 
current understanding of binary and stellar evolution, 
arguing strongly that the maximum neutron star mass is 
below 2.3-2.4 M 0 . Advanced LIGO will be able to place 
more stringent constraints on this merger rate, and this 
measured merger rate will provide more firm constraints 
on the maximum neutron star mass. 

Unfortunately, constraining the NS-NS merger rate 
and its comparison with short GRB rates may not di¬ 
rectly constrain the equation of state. BH-NS merg¬ 
ers can also produce short-duration GRBs and, in a 
small subset of binary population synthesis scenarios, 
BH-NS merger rate is on par with the NS-NS merger 
rate, e.g.(Dominik et al. 2012). Although it is possible 
to construct a scenario where BH-NS mergers dominate 
the rate of NS-NS mergers, recall from Section 2.1 that 
it is likely that most of BH-NS systems do not p roduce 
disks of sufficient size to make GRBs (Rosswog 2015). 
If the connection between the short GRBs and BH-1NS 
mergers could be established, it will place strong con¬ 
straints on binary population synthesis models (in most 
scenarios, the NS-NS merger rate is an o rder o f magni- 
tude higher than the BH-NS merger rate; [Dominik et al.[ 
(2012). It will also point to the EOS with maximum NS 
mass over 2.3-2.4M 0 as such EOS allows to eliminate 
the formation of short GRBs from NS-NS mergers via 
most accepted (BHAD) engine model. 

Alternatively, it is possible that the baryon loading 
from neutrino-driven winds from neutron stars does not 
preclude the possibility that cores that remain neutron 
stars can produce GRBs. Late-time emission from short- 
duration GRBs either argues for some extended power 
source (e.g. late- time accr etion, magnetar) or jet inter¬ 
act ion (Holcomb et al. 2014). The magnetar power source 
requires a long-lived neutron star and has been shown, 


within the free pa rameters of the magnetar model (Rowl- 


inson et al. 12014 1, to fit the late-time emission well. .Pro¬ 


ponents ot the magnetar model argue that there must be 
some way to avoid the burying of the magnetic field and 
the baryon loading, e.g. a high pressure, magnetic or ra¬ 
diation dominated, region along the axis prevents baryon 
loading and allows the production of a relativistic jet. 

To rule out these alternative solutions, we either need 
to refine our theoretical models to prove the differences in 
the current interpretations or we need additional obser¬ 
vational diagnostics to distinguish the above picture. For 
example, it might be possible to distinguish BH from NS 


formation using the gravitational wave signal (Chatzi- 


ioannou et al. 20151, but this will likely require a very 
strong signal from a nearby merger. If a stable neutron 
star with strong magnetic fields is formed, we would ex¬ 
pect continued pulsar-like activity, e.g. a soft gamma-ray 
repeater. However, if we only have LIGO localizations, 


it will be difficult to pr ove t he associati on of any activity 


with the merger, e.g.(Kelley et al. 2010). The ejecta from 
these models could allow us to distinguish these models. 
The optical/infra-red light produced by the decay of ra¬ 
dioactive elements can be used to distinguish between NS 
and BH fates. The masses of the dynamic ejecta that are 
shown in figurejo]) do not vary dramatically in our models 
and, without alteration by the core, the extreme neutron- 
richness of these ejecta leads to a composition that al- 


large opacities (Kasen et al. 

2013) producing transients 

peaking in the near-in 

rared (or infrared) (Roberts et al. 

2011; Barnes & Kasen 

2013i|Tanaka & Hotokezaka|2013; 

Grossman et ai.]|2014 

Fontes et al.||2015|. However, the 


produced by a combination of this dynamic ejecta, neu¬ 
trino driven winds from the core (if the core does not 
collapse to a black hole) and ejecta from the disk accre¬ 
tion. Due to the longer exposure to the intense neutrino 
background, the neutron fraction is reduced, altering the 
composition of these ejecta and the revised opacities may 
allow bri ghter transien ts in the near-infrared and opti¬ 


cal) Perego et al. 2014). Since this ejec ta is much larger 


for IMS systems than B H disk systems)Fernandez et al. 


2015 Just et al. 20151) , we may be able to distinguish 


these two systems by the wavelength of their electromag¬ 
netic emission. If so, we can use observations of neutron 
star mergers to place constraints on dense matter equa¬ 
tions of state. 



ejecta mass [M 0 ] 

Figure 6. The number of systems as function of ejecta mass in 
the tidal tails from these neutron star mergers for our increased 
NS models for both solar and l/10th solar metallicity. Note that 
althought the ejecta masses range from below 0.01 Mq to above 
0.03 Mq, 75% of the mergers produce ejecta with masses of 0.025 =t 
0.003 M©. 
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Table 1 

Merger Properties * 1 


M NS i 

(Me) 

MnS 2 

(M e ) 

M b 

iV1 core 

(Me) 

Q b 

►~>core 

(k B /nuc) 

I b 

Jcore 

( 10 49 gcm 2 s -1 ) 

M disk b 

(Me) 

]\ \ . b 

iVA ejected 

(Me) 

T equator C 

(km) 

^pole 

(km) 

1.0 

1.0 

1.51 

0.91 

1.43 

0.48 

0.0084 

21.3 

10.9 

1.2 

1.0 

1.58 

0.93 

1.12 

0.59 

0.029 

20.9 

12.9 

1.4 

1.0 

1.74 

0.94 

1.22 

0.63 

0.028 

21.9 

13.5 

1.6 

1.0 

1.86 

1.07 

1.20 

0.71 

0.030 

18.2 

13.6 

1.8 

1.0 

2.10 

1.06 

2.27 

0.67 

0.032 

71.1 

14.2 d 

2.0 

1.0 

2.14 

1.28 

1.41 

0.83 

0.032 

175.3 

14.4 d 

1.2 

1.2 

1.87 

0.94 

1.91 

0.52 

0.010 

22.4 

12.6 

1.4 

1.2 

2.01 

0.95 

1.92 

0.57 

0.022 

24.0 

13.4 

1.6 

1.2 

2.10 

0.97 

1.88 

0.67 

0.034 

23.2 

14.0 

1.8 

1.2 

2.19 

1.19 

1.87 

0.77 

0.036 

21.9 

14.9 

2.0 

1.2 

2.31 

1.43 

1.89 

0.85 

0.037 

73.5 

14.6 d 

1.3 

1.4 

2.08 

1.01 

2.22 

0.61 

0.016 

21.6 

12.7 

1.3 

1.4® 

2.02 

0.69 

1.86 

0.67 

0.016 

20.6 

12.9 

1.4 

1.4 

2.29 

0.96 

3.11 

0.50 

0.012 

25.2 

11.6 

1.6 

1.4 

2.33 

0.95 

2.41 

0.65 

0.021 

25.4 

13.9 

1.8 

1.4 

2.40 

1.00 

2.15 

0.76 

0.038 

25.2 

14.3 

2.0 

1.4 

2.48 

1.45 

2.15 

0.88 

0.041 

25.4 

14.7 

1.6 

1.6 

2.55 

1.03 

3.07 

0.64 

0.013 

23.7 

13.6 

1.8 

1.6 

2.71 

1.08 

3.36 

0.67 

0.019 

25.2 

13.8 

2.0 

1.6 

2.70 

1.36 

2.76 

0.86 

0.042 

24.5 

14.6 

1.8 

1.8 

2.87 

1.08 

3.84 

0.71 

0.015 

25.9 

13.6 

2.0 

1.8 

3.04 

1.28 

4.11 

0.74 

0.019 

26.4 

14.2 

2.0 

2.0 

3.31 

1.19 

5.37 

0.67 

0.017 

14.1 

13.3 


a All models from ( [Korobkin et al.|2012[ ) with the simulation physics from ( [Ros swog 2013|). 

b The core is defined by the post-merger material whose density is above 10 l 4 gcm —a at the end of the 

calculation. 


c ^*pole 5 ^equator correspond to the radius of the compact core along the rotation axis and equator re¬ 
spectively assuming the core reaches a uniformly-rotating, equilibrium state using the FSU2.1 equation 
of state. 

d For the most part, the angular momenta in these proto-neutron star is just below most limits for 
secular instabilities. For these 3 systems, this is not the case. It is possible that these instabilities can 
quickly shed angular momentum. 

e This model assumes co-rotation for the initial conditions. 
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Table 2 

Remnant Fate 


MnSI 

M N S2 

Fate (eq[2|| a 

Fate (eq[3]l a 

Fate (eq[3|| a 

Fate (eq[3l a 

Fate (eq[3]) a 

Fate (eq|3| a 

Fate (eq|3|l a 

(M 0 ) 

(M 0 ) 

FSU2.r 

M 0 = -h 5 

A4q = zTz 

M 0 = ±i 

M 0 = Ti 

M 0 = 

Mq = z77 


1.0 

1.0 

NS 

NS 

NS 

NS 

NS 


NS 

NS 

1.2 

1.0 

-BHspin 

BHacc : 210 

NS 

NS 

NS 


NS 

NS 

1.4 

1.0 

BH acc : 180 

BHacc : 40 

BHacc : 220 

BHacc : l-2s 

NS 


NS 

NS 

1.6 

1.0 

BH acc : 90 

BHacc : 90 

BHacc : 90 

BHacc : 100 

BHacc : 

240 

BHacc : 1.8s 

NS 

1.8 

1.0 

BH a cc : 60 

BH 

BHacc : 20 

BHacc : 90 

BHacc : 

120 

BHacc : 360 

^H s pin 

1.8 

1.0 b 

BH 

BH 

BH 

BH 

90 


200 

BH s pi n 

2.0 

1.0 

BHacc : 50 

BH 

BH 

BHacc : 50 

BHacc 

50 

BHacc : 50 (BH) 

BHacc : 170 

2.0 

1.0 b 

BH 

BH 

BH 

BH 

BH 


BH 

90 

1.2 

1.2 

BHacc : 500 

BHacc : 50 

BHacc : 450 

BHacc 3s 

NS 


NS 

NS 

1.4 

1.2 

BH a cc : 50 

BH 

BHacc : 50 

BHacc : 280 

BHacc : 

480 

BH sp i n 

NS 

1.6 

1.2 

BH a cc : 10 

BH 

BHacc : 20 

BHacc : 30 

BHacc 

70 

BHacc : 200 

BHspin 

1.8 

1.2 

BHacc : 10 

BH 

BHacc : 20 

BHacc : 30 

BHacc 

30 

BHacc : 30 

BHacc : 200 

2.0 

1.2 

BH 

BH 

BH 

BH 

BH 


BHacc : 60 

BHacc : 60 

2.0 

1.2 b 

BH 

BH 

BH 

BH 

BH 


BH 

BHacc : 30 

1.3 

1.4 

BH a cc : 60 

BHacc : 60 

BHacc : 60 

BHacc : 70 

BHacc : 

300 

BHacc 3s 

NS 

1.3 

1.4 C 

BHacc • ~ 3s 

BHacc • ^ 3s 

BHacc • ^ 3s 

BHacc • ^ 3s 

BHacc : 

~ 3s 

BHacc ■ ~ 3s 

NS 

1.4 

1.4 

BHacc : 60 

BH 

BHacc : 60 

BHacc : 60 

BHacc 

60 

BHacc : 490 

BHspin 

1.6 

1.4 

BH 

BH 

BH 

BH 

BHacc 

30 

BHacc : 60 

BHacc : 700 

1.8 

1.4 

BH 

BH 

BH 

BH 

BHacc 

20 

BHacc : 40 

BHacc : 80 

2.0 

1.4 

BH 

BH 

BH 

BH 

BH 


BHacc : 40 

BHacc : 40 

1.6 

1.6 

BH 

BH 

BH 

BH 

BH 


BHacc : 30 

BHacc : 420 

1.8 

1.6 

BH 

BH 

BH 

BH 

BH 


BHacc : 40 

BHacc : 130 

2.0 

1.6 

BH 

BH 

BH 

BH 

BH 


BH 

BHacc : 50 

1.8 

1.8 

BH 

BH 

BH 

BH 

BH 


BH 

BHacc : 70 

2.0 

1.8 

BH 

BH 

BH 

BH 

BH 


BH 

BHacc : 30 

2.0 

2.0 

BH 

BH 

BH 

BH 

BH 


BH 

BHacc : 40 


a The final fates are NS, the merged system forms a NS AD that, even after the accretion of a disk, remains a NS; BH a cc, merged 
system initially forms a NS AD that, after some accretion, collapses to a BH forming a BHAD (the accretion time is given in 
ms unless otherwise specified); BH sp i n , merged system initially forms a NS AD that only collapses after the disk accretion and 
the loss of angular momentum; BH, the merged system collapses immediately to a BH, forming a BHAD; BH sp i n (the merged 
system remains a NS AD until the disk accretion finishes after which spin down is able to decreases the maximum neutron star 
mass and allow it to collapse to a black hole.) 

b For these models, we consider the fate of the neutron star with a reduced rotation assuming the secular instabilities quickly 
reduce the angular momentum. Any changes are noted in parantheses. Note, however, that typical growth timescales for these 
systems are estimated to be above 100 ms and it is more likely that the spin down occurs after the accretion phase. 
c This model assumes co-rotation for the initial conditions. 
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Table 3 

Merger Fate Distributions a 


Equation of State b 

BHAD GRB a 

BH a 

BHacc, tacc < 100ms a 

BH a cc, t a cc ^ 100mS a 

Dtr a 

^-n-spin 

NS a 

FSU2.1 C 

0.8739 

(0.8715) 

0.0106 

(0.0100) 

0.8633 

(0.8615) 

0.1211 

(0.1259) 

0.0051 (0.0026) 

0.0000 

(0.0000) 

Steiner2.0 

0.9988 

(1.0000) 

0.8730 

(0.8707) 

0.1258 

(0.1293) 

0.0012 

(0.0000) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner 2.2 

0.8723 

(0.8707) 

0.0090 

(0.0091) 

0.8633 

(0.8615) 

0.1209 

(0.1259) 

0.0067 (0.0035) 

0.0000 

(0.0000) 

Steiner2.3 

0.0510 

(0.0786) 

0.0104 

(0.0095) 

0.0406 

(0.0690) 

0.8734 

(0.9010) 

0.0532 (0.0139) 

0.0224 

(0.0065) 

Steiner2.4 

0.0189 

(0.0113) 

0.0066 

(0.0087) 

0.0123 

(0.0026) 

0.1514 

(0.0924) 

0.5543 (0.4952) 

0.2754 

(0.4010) 

Steiner 2.5 

0.0104 

(0.0095) 

0.0032 

(0.0048) 

0.0072 

(0.0048) 

0.0104 

(0.0013) 

0.0671 (0.0694) 

0.9120 

(0.9197) 

Steiner2.7 

0.0034 

(0.0048) 

0.0000 

(0.0000) 

0.0034 

(0.0048) 

0.0057 

(0.0048) 

0.0072 (0.0009) 

0.9838 

(0.9896) 

FSU2.1 d 

0.8704 

(0.8662) 

0.0065 

(0.0040) 

0.8639 

(0.8622) 

0.1248 

(0.1323) 

0.0048 (0.0015) 

0.0000 

(0.0000) 

Steiner2.0 

0.9982 

(1.0000) 

0.8690 

(0.8642) 

0.1291 

(0.1358) 

0.0018 

(0.0000) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner 2.2 

0.8694 

(0.8662) 

0.0055 

(0.0040) 

0.8639 

(0.8622) 

0.1253 

(0.1319) 

0.0053 (0.0020) 

0.0000 

(0.0000) 

Steiner2.3 

0.0414 

(0.0652) 

0.0063 

(0.0040) 

0.0351 

(0.0612) 

0.8839 

(0.9210) 

0.0568 (0.0094) 

0.0179 

(0.0044) 

Steiner2.4 

0.0134 

(0.0074) 

0.0040 

(0.0040) 

0.0094 

(0.0035) 

0.1587 

(0.0785) 

0.5610 (0.5042) 

0.2669 

(0.4099) 

Steiner 2.5 

0.0065 

(0.0040) 

0.0012 

(0.0010) 

0.0053 

(0.0030) 

0.0087 

(0.0015) 

0.0631 (0.0612) 

0.9217 

(0.9333) 

Steiner2.7 

0.0015 

(0.0005) 

0.0000 

(0.0000) 

0.0015 

(0.0005) 

0.0042 

(0.0035) 

0.0052 (0.0010) 

0.9891 

(0.9951) 

FSU2.1 e 

0.8981 

(0.8724) 

0.0159 

(0.0100) 

0.8821 

(0.8624) 

0.1019 

(0.1276) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner2.0 

1.0000 

(1.0000) 

0.8851 

(0.8355) 

0.1149 

(0.1645) 

0.0000 

(0.0000) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner 2.2 

0.8964 

(0.8724) 

0.0144 

(0.0100) 

0.8820 

(0.8624) 

0.1036 

(0.1276) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner2.3 

0.8923 

(0.8724) 

0.0149 

(0.0100) 

0.8774 

(0.8624) 

0.1077 

(0.1276) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner2.4 

0.0566 

(0.0799) 

0.0107 

(0.0087) 

0.0458 

(0.0712) 

0.8726 

(0.8655) 

0.0708 (0.0547) 

0.0000 

(0.0000) 

Steiner 2.5 

0.0159 

(0.0100) 

0.0055 

(0.0074) 

0.0104 

(0.0026) 

0.0523 

(0.0438) 

0.8182 (0.7947) 

0.1136 

(0.1515) 

Steiner2.7 

0.0052 

(0.0078) 

0.0000 

(0.0000) 

0.0052 

(0.0078) 

0.0086 

(0.0022) 

0.0164 (0.0039) 

0.9698 

(0.9861) 

FSU2.1 f 

0.8961 

(0.8667) 

0.0107 

(0.0044) 

0.8854 

(0.8622) 

0.1039 

(0.1333) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner2.0 

1.0000 

(1.0000) 

0.8839 

(0.8326) 

0.1161 

(0.1674) 

0.0000 

(0.0000) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner 2.2 

0.8938 

(0.8667) 

0.0085 

(0.0044) 

0.8853 

(0.8622) 

0.1062 

(0.1333) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner2.3 

0.8886 

(0.8662) 

0.0100 

(0.0040) 

0.8786 

(0.8622) 

0.1114 

(0.1338) 

0.0000 (0.0000) 

0.0000 

(0.0000) 

Steiner2.4 

0.0459 

(0.0652) 

0.0062 

(0.0044) 

0.0398 

(0.0607) 

0.8806 

(0.8785) 

0.0735 (0.0563) 

0.0000 

(0.0000) 

Steiner 2.5 

0.0107 

(0.0044) 

0.0038 

(0.0030) 

0.0068 

(0.0015) 

0.0478 

(0.0360) 

0.8240 (0.7985) 

0.1176 

(0.1610) 

Steiner2.7 

0.0023 

(0.0020) 

0.0000 

(0.0000) 

0.0023 

(0.0020) 

0.0058 

(0.0025) 

0.0167 (0.0049) 

0.9751 

(0.9906) 


a We include 4 suites of simulations corresponding t wo different initial conditions: standard/old StarTrak initial conditions(|Belczynski et al.| 
|2008|) and the new initial conditions (|Sana et al.|2012|) and two different initial compact remnant mass distributions:. For each result, we include 
two values, one where all systems are included and one where it is assumed that the binary is destroyed if the system goes through a common 
envelope in the Hertzprung Gap. For the remnant fates, we consider those systems that collapse immediately (within 3 ms of the merger) to a 
BH (BH), systems that collapse to a BH within 100 ms (BH acc , tacc < 100ms), those that collapse to a BH through accretion taking more than 
100ms (BHacc, t acc > 100ms), those that collapse only after spin-down (BH sp i n ), and those that remain neutron stars (NS). The first column 
is the sum of the second and third columns (BH+BH acc , tacc < 100ms) and represents those systems that can form the canonical black hole 
accretion disk GRB without significant baryonic contamination. 

b The equations of state correspond to the FSU2.1 equation of state with a maximum non-rotating neutron star mass of 2.1 Mq and a range of 
equations of state from our parameterized model (SteinerX) where “X” denotes the maximum non-rotating neutron star mass for that equation 
of state (in Mq). 

c This suite studies the old initial conditions with the standard NS mass distribution. 
d This suite studies the new conditions with the standard NS mass distribution. 
e This suite studies the old initial conditions with the enhanced (+0.1 Mq) NS mass distribution. 
f This suite studies the new new conditions with the enhanced (+0.1 Mq) NS mass distribution. 













